Dynamic generation of entangling wave packets 
in XY spin system with decaying long range couplings 
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We study the dynamic generation of spin entanglement between two distant sites in an XY model 
with l/r'^-decay long range couplings. Due to the linear dispersion relation e(fc) ~ of magnons 
in such model, we show that a well-located spin state can be dynamically split into two moving 
entangled local wave packets without changing their shapes. Interestingly, when such two wave 
packets meet at the diametrically opposite site after the fast period r = vr, the initial well-located 
state can be recurrent completely. Numerical calculation is performed to confirm the analytical 
result even the ring system of sizes up to several thousands are considered. The truncation 
approximation for the coupling strengths is also studied. Numerical simulation shows that the 
above conclusions still hold even the range of the coupling strength is truncated at a relative shorter 
scale comparing to the size of the spin system. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Ud, 75.10.Jm 



INTRODUCTION 

In quantum information processing, it is also crucial 
to generate entangled qubits, which can be used to per- 
fectly transfer a quantum state over long distance. For 
optical system this task has been completed long time 
ago, but for a solid state systems it remains a great chal- 
lenge both in experiment and theoretical setup to create 
quantum entanglement by a solid state device. Recently 
many proposals to entangle distant spins have beenpro- 
posed based on various physical mechanisms [ll S S H ■ 
Most protocols for accomplishing quantum state transfer 
in a spin array base on the fixed int er-q ubit couplings 

liilliiEIIlIIlIIlIIlElIIlIIlIIi The Sim- 
plest coupled spin system with uniform nearest neighbor 
(NN) couphng has been studied by the pioneer work 0. 
The paradigm to generate maximal entanglement and to 
perform perfect qubit-state transmission over arbitrary 
distance is the protocols using the pre-engineered inho- 
mogeneous NN couplings Among them, the prac- 

tical scheme realizes the transmission of Gaussian wave 
packet as flying qubit via the spin system with uniform 
NN couplings |3, 0|. The advantages of such scheme 
rest with its fast transfer, which means that the period 
of transfer time is proportional to the distance, and the 
higher fidelity for longer distance. We also notice that 
there is a current interest in stud y in g systems with long- 
range inter-qubit interactions [2iJ, |2l|, |22| . 

In this paper, we revisit the issue of entanglement 
generation and quantum state transfer between two dis- 
tant qubits in a qubit array with long-range inter-qubit 
interactions. An alternative way to construct a per- 
fect medium for quantum state transfer may require the 
long-range interactions beyond NN interactions, but they 
should be required to decay rapidly in order to avoid the 
direct connecting coupling between two distant qubits 
that trivially causes quantum entanglement. Now we 



FIG. 1: (color online) Schematic illustration for the ideal dis- 
persion relation e{k) ~ |fc| of a system which is shown to be a 
perfect entangler. Such kmd of system can be realized by the 
XY model with pre-engineered long-range couplings. 



propose a novel protocol based on a pre-engineered XY- 
model with long-range 1/r^-decay interactions. It is 
found that such model has the same function as that 
of the modulated NN couphng spin model 0, offering 
significant advantages over other protocols in the tasks 
perfectly transferring quantum state and generating en- 
tanglement between two sites over longer distance. 

This paper is organized as follows, in Sec. II, we con- 
sider a model with linear dispersion relation and show 
that it ensures the perfect state transfer and creation of 
entanglement between spatial separated qubits. In Sec. 
Ill we propose a spin model with 1/r^-decay interactions 
that process the desired dispersion relation. In Sec. IV 
we perform numerical calculation to confirm some of the 
obtained analytical results and the validity of the trun- 
cation approximation. 
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FIG. 2: (color online) Schematic illustration of the generation 
process of two distant entangled qubits. At t = the initial 
state is a single-spin flip on the saturated ferromagnet. Att > 
the initial 6 -pulse magnon is separated into two components 
which have maximal entanglement. At t = t, the recurrence 
time, the two sub-wave packets meet together and the initial 
state recurrent. 



FORMALISM 
Pre-engineered model with linear dispersion relation 

Usually, we regard a 1/2 spin as a qubit, and a coupled 
spin system as a qubit array respectively. The simple 
coupled qubit array is usually described by the spin- 1/2 
XY model. Our proposal makes the quantum spin array 
behave as a perfect spin- networks, in which a spin flip 
at any site (or the superposition of local single-spin-flip 
states) can evolve into two entangled, local wave pack- 
ets. Consider the prescribed spin-1/2 XY-mode\ on ring 
system with N sites. The Hamiltonian reads 



(1) 



(2) 



the single excitation (magnon) subspace, the basis states 
are denoted by the single-site flipped state (or ^-pulse) 
|i) = |0), i e [1,-/V]. Therefore, if state |i) can evolve 
to state \j) after a period of time r, \j) = exp{—iHT) 
we have 



N 



N 



(«IT),+/3U),)niT); 
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"("IT),:+/3U)jnit)M 
(3) 

i.e., a qubit state a|t) + can be transferred from 
i to j. Furthermore, any single-magnon state can be 
expressed as 



'A,H) 



and 



k 



(4) 



(5) 



respectively in spatial and momentum {\k)} spaces. 
Here, we have used the spin-wave operator 



(6) 



with discrete momentum k = 2TTn/N, n e [—N/2, N/2 — 
1]. 

We begin with assumption that there exists an opti- 
mal distribution of J^, which ensures the single-magnon 
spectrum possessing a linear dispersion relation, i.e. 
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k 

N 
2^ 



(7) 



|fc| 



as illustrated in Fig. 1. In the following, we first show 
that such kind of systems can perform perfect state trans- 
fer and long range entanglement generation, and then we 
can provide a practical example that satisfies this linear 
dispersion relation. 



where 5f , Sf and Sf are Pauli matrices for the site «th, 
and sf = Sf ± iSf; Jr is the coupling strength for the 
two spins separated by the distance r. Since the Hamil- 
tonian conserves spin, i.e., [^,- S^, H] = 0, the dynamics 
can be reduced to that in some invariant subspaces. Thus 
we can only concentrate on the single excitation subspace 
hereafter. If a single-site flipped state can be correctly 
transferred, a qubit state should also be transferred corre- 
spondingly because the saturate ferromagnetic state with 
all spins up state |0) = Hi^i is an eigenstate of the 
Hamiltonian. Actually, in the single-spin-flip subspace or 



Time evolution of wave packets in the linear 
dispersion regime 

In order to investigate the dynamics of generating en- 
tanglement in the spin system with linear dispersion 
relation mentioned above, we concentrate on the case 
where the initial state is single-site flipped state S~ |0), 
i £ [1,^]- Intuitively, a well locahzed state has equal 
probability amplitudes with respect to momentum eigen- 
states I ± I A: I), due to the inverse Fourier transformation 
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integer. In the following, we only consider the states 
at discrete instant t, which ensures a to be integer, i.e., 
Nt/2TT = [Nt/2Tr] (the integer part of Nt/2TT). In terms 
of 



N 
2^' 



(15) 



Thus it should be split into two classes of waves 
with ± |fc| driven by the Hamiltonian with linear single- 
niagnon dispersion relation 0. If the superposition of 
the two classes of waves are still well-localized (or form 
two wave packets) in real space, it will lead to the entan- 
glement of the two wave packets. 

Generally we first consider a wave packet dynamics 
with an initial state 

mNA,0)) = 5^^ |0) = |0) , (9) 

which is a single flip at site Na- The time evolution 
starting with this state can be calculated as 



(8) the wave packets H12|l can be re-written as 

\MNa, t)) = ^ ± 1]} \j) . (16) 



Obviously, Eq. H14|l shows that state \4>±{NA,t)) is 
shape- invariant with velocity v± = ±iV/27r. Now we 
show the locality of these states. The projections 



{j \<t>±{NA,t)) 



odd l±{j,t) 



0, evenl±{t)^0 



(17) 



of \<p±{NA,t)) onto the localized state \ j) shows that the 
shape-invariant states are well-localized. And each wave 
packet has the probability amplitude 



-iHt 



^NaM 



\d,.{NA.t)) + \^+{NA,t)) (10) 
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n=0 
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where 



mNA,t)) = £\^''(^-^^'Uk\k) 



(11) 



denote the two wave packets with explicit forms in the 
real space 



which indicates that the initial state © is splited into 
two wave packets completely. On the other hand, the 
dispersion |(7J) requires the translational invariance of the 
Hamiltonian since the momentum is the conserved quan- 
tity, which leads to the periodicity of the states |(/)±) 



\MNA,t)) 



2^ Jo Vn Y 
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(12) r = 



where m± — 1,2,- ••, are integers. Obviously, when 
T — {171+ + m_)7r, the initial state |i/)(iV^,0)) = 5*^^ |0) 
evolves to 



We will show that such two wave packets are localized 
and nonspreading, with velocity v± = ±iV/27r, respec- 
tively. Actually, from Eq. H12(l one have 



|0) 



(20) 



\(t>±{NA,t + to)) 



N 

c^±{NA±—t,to)) (13) 



or 



-iHt 



\^±iNA,to))^T.,\^±{NA,to)) (14) 



where Ta acts as translational operator for the states 
|(/)±(Z,0): Ta\(l)±{l,t)) = \(j)±{l±a,t)). In general, the 
spatial coordinate is treated as discrete variable while 
the temporal coordinate is continuous. In order to make 
Ta to be operable, the translational spacing a should be 



which is just the recurrence of \'iI){Na, 0)) on the positions 
Na — N{m+ — m_)/2 at instants (m+ -I- m-)'!T. The 
physics of this phenomenon can be understood as the 
interference of two wave packets H12|l . It also accords with 
the prediction for the system with spectrum-symmetry 
matching condition (SSMC) El [13. 

In the single magnon invariant subspace, the single 
flip states {\j)} at site j constitute the complete ba- 
sis. Then the above conclusion can be applied to any 
states in this subspace. Consider an arbitrary initial state 
1^(0)) = \j) ^ A, |V'(j, 0)), which is a coherent 

superposition of single flip states. At time t, it evolves to 



\m) 



$(0)) = |$+(i)) + |$_(t)> 



(21) 
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where 



(22) 



where i?=even (odd) represents the parity of the state 
under the reflection. The concurrence between two spins 
at Na + j and Na - j for the state \^{t)) is 



represent two invariant-shape states. Accordingly, at in- 
stants r = (m+ -|- m-)'!T, the final state is a translation 

of \m) 



(m_j_ — m_ ) 



\m) 



(23) 



Furthermore, during the period of time t =/= t, wave pack- 
ets |$±(i)) are still well localized in space if the initial 
state |$(0)) is local. Then the revival of a wave packet 
can be used to implement perfect quantum state transfer. 



Entanglement of two separated spins 

Now we turn our attention on the entanglement of the 
two separated spins induced by the wave packets p2|l . 
The reduced density matrix of a state |<&(t)) for two spins 
located at sites i and j |2^ has the form 



l<3 

w 









V 



Z,j w 





\ 



^J 

-^J 



(24) 



with respect to the standard basis vectors I tt), |Ti)i I it); 
and lii). Here, the matrix elements are 



Wi 



m)\{S^S;+SrSt)mt))- (25) 



Correspondingly, the concurrence of two spins located at 
sites i and j for a state \^{t)) can be calculated by 



Cij — max < 0, 2 



(26) 



Since the state we concern is in the invariant subspace 
with Sz = N/2 — 1, we have v^^ = and then the con- 
currence reduces to 



c,, = \{m\{sts- + srsi)\m)\ 



(27) 



Consider an initial state having reflectional symmetry 
with respect to a point Na, which has the form 



(29) 



Obviously, the total concurrence ^jyoC{j, t) = 1 due to 
the normality of the wave function |$(i)). Usually, the 
total concurrence is regarded as the concurrence of two 
wave packets. Taking the single flip at A''^ in the form 
Q to be the initial state as an example, straightforward 
calculation shows that non-zero concurrences arc 



(2n - 1)2^2 



(30) 



For this to occur, one requires that to = 2Trj/N, tn = 
27r[j + (2n - l)]/iV, n = 1,2, .... It indicates that the 
concurrence with magnitudes (|3()|l for two spins at Na+J 
and Na — j can be generated at the moment t^, ti, t2, 



MICROSCOPIC LONG-RANGE INTERACTION 
MODEL 



In this section, we consider the possibility to realize 
the above scheme in a system which possesses the linear 
dispersion relation based on an one-dimensional arrange- 
ment of spins (qubits) coupled by long-range interactions. 
Actually, using the identity 



N 



E 

r—odd 



AN r7r(2£) 
^T-^ cos ■ 



N 



(31) 



we have 



N 
2^ 



N/2-l 



Jo + ^ 2JrCos{kr), 



(32) 



r—odd 



where 



\m)-fNA0)\NA)+J2fjmNA+j) + {-lf\NA-j)], 4' rV ^''''> 

(28) Then the Hamiltonian matches the case {7} in single- 
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magnon invariant subspace can be rewritten as 

N/2-l 
r—odd 

ho - TV Jo, 

hr = J2{S+Sr^, + h.c.). (34) 

i 

Notice that, although such model involves the long range 
interaction, Jr decays rapidly as r increases. So we 
call this as "microscopic long-range interactions" . The 
above analytical analysis shows that such microscopic 
long-range interaction can lead to nontrivial long-range 
entanglement and QST. 



NUMERICAL RESULTS AND TRUNCATION 
APPROXIMATION 

In the previous sections, it is found that the micro- 
scopic long-range interaction can lead to nontrivial long- 
range entanglement and QST for large N limit system. 
In this section, numerical simulations are performed for 
finite systems to illustrate the result we obtained above 
and investigate the application of it. The numerical ex- 
act diagonalization method is employed to calculate the 
time evolution of the single flip state (jH)) and the Gaus- 
sian wavepacket in the following form 



\i^GiNA,a)) = ^e-^(^--)^5r |0) (35) 

in the finite N system. Here, il is the normalized fac- 
tor and Na, a determine the center and the shape of 
wavepacket. In Fig. 3 and 4, the time evolution of a 
single flip state ^ and a Gaussian wave packet (|35(l of 
a = 0.1 in a iV = 100 ring system are plotted. It shows 
that the local initial states split into two wave packets in 
both cases and the wave packets keep their shapes with- 
out spreading. Similarly, numerical simulations are also 
employed for the open chain systems. The similar con- 
clusions are also obtained for open chain system which 
will be discussed in the following. 

A natural question is that what role the long-range 
coupling plays for the peculiar behavior of the propaga- 
tion of the local wave in such a system. As shown above, 
the long-range hopping coupling decays rapidly. Then 
if the size of system is large enough, it is believed that 
the too long range couphng can be neglected. Thus, the 
"long range coupling" in the form (p?!^ can be regarded as 
the relative local, or "microscopic long-range" coupling. 
In order to investigate this problem, or the boundary be- 
tween the so called "long range" and "microscopic long- 
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FIG. 3: (color online) Time evolution of a 5-pulse in a N = 
100 ring system obtained by numerical simulation. It shows 
that the two sub-waves are local and keep the shape. 
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FIG. 4: (color online)Time evolution of a Gaussian wave 
packet of a — 0.1 in a N — 100 ring system obtained by 
numerical simulation. It shows that the two sub-waves are 
local and keep the shape. 

range" couplings, we consider the truncated Hamiltonian 
H = ho+ JrK, (36) 

r—odd, <ro 

where tq is the truncation distance. From the above 
analysis, the initial state \tp{NA,0)) should recurrence 
at the positions Na at instant 27r, if the system is per- 
fect for the quantum state transfer. The autocorrelation 
|^(i)| = \{i{}{Na,0) \tp{NA,t))\ is appropriate quantity 
to investigate the role that rg plays. During the pe- 
riod of time ~ [0,37r], the maxima of autocorrelations 
l^maxl — ina.x{\A{t)\} of the state © in the systems 
with N = 500, 1000, 1500, 2000, and ro = 10, 20, 100 
are calculated numerically by the exact diagonalization 
method. In Fig. 5, the dependence of the truncation 
distance rg on the quantity |j4i„ax| is plotted. It shows 
that |Amax| approaches to 1 when rg is around 90, which 
is called the critical truncation distance or the boundary 
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FIG. 5: (color online) The maximal autocorrelation functions 
of the initial 5 -pulse state in the N-site systems with the trun- 
cation distance ro. It demonstrates that critical ro, at which 
l^maxl start to approach to 1, do not depend on N strongly. 



FIG. 6: (color online) The time- dependent concurrences of 
two spins separated by the distance 21 = 800 in the chain 
system with N = 1000, ro = 10,15,20,30 and 500. The re- 
sults for ro = 500 is in agreement with the analytical analysis. 
The results for different truncation approximations show that 
range of LR coupling can be taken in a small scale due to the 
I /r^ -decay of coupling constants. 



between the "long range" and "microscopic long-range" 
couplings, for the systems with different N. It indicates 
that in the case of iV ^ tq, the wave packets still travel 
without spreading. By making use of this observation, 
we find that although the interactions between spins are 
"long range", the 1/r^-distribution of coupling strength 
allow us to limit the maximal interaction range while 
minimizing the degradation of the quantum coherence 
obtained for ideal model. Then we have the conclusion 
that the l/r'^-decay coupling can be regarded as local 
coupling, or microscopic long-range coupling. In other 
word, the robust long-range entanglement between two 
distant qubits is not due to the direct long range cou- 
pling interaction between them. 

From the above observation, we find that, for large 
size system the long-range coupling (ro > 90) can be ne- 
glected. Then in thermodynamic limit, a ring system is 
equivalent to a more practical system, open chain sys- 
tem. To demonstrate this, the numerical simulation is 
employed to investigate the concurrence C{l,rQ,t) be- 
tween two far separated sites N/2 zL I {I ^ -^/2) for the 
TV-site system with different truncation rp. In Fig. 6, 
plots of C{1, ro,t) for the systems with N = 1000, / = 400 
and different tq are presented. It shows that for an open 
chain system, long range entanglement between two dis- 
tant qubits can be achieved via the time evolution of a 
single flip state. The maximal entanglement created by 
such system is 0.5, as measured by the two-point concur- 
rence. 



SUMMARY 

In summary, the system with long-range coupling is in- 
vestigated analytically and numerically. It is found that 
the 1 /r^-decay long-range coupling model can exhibit ap- 
proximately linear dispersion e \k\. The dynamics of 
such model possesses a novel feather that an initial local 
wave packet can be separated into two entangled local 
wave packets. Furthermore, during the traveling period 
each wave packets can keep their shapes without spread- 
ing. Numerical simulation indicates that there exists a 
critical truncation distance tq, which limits the range of 
the interaction but not affects the generation of entangle- 
ment between two distant qubits in the distance / ^ tq. 
This model open up the possibility to realize the solid- 
state based entangler for creating two entangled but spa- 
tially separated qubits. 
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